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We present a systematic study that clarifies validity and limitation of current hybrid function- 
als in density functional theory for structural and electronic properties of various semiconductors 
and insulators. The three hybrid functionals, PBEO by Perdew, Ernzerhof, and Becke, HSE by 
Heyd, Sucseria, and Ernzerhof, and a long-range corrected (LC) functional, are implemented in 
a well-established plane-wave-basis-set scheme combined with norm-conserving pseudopotentials, 
thus enabling us to assess applicability of each functional on equal footing to the properties of the 
materials. The materials we have examined in this paper range from covalent to ionic materials as 
well as a rare-gas solid whose energy gaps determined by experiments are in the range of 0.6 eV - 
14.2 eV: i.e., Ge, Si, BaTi0 3 , /3-GaN, diamond, MgO, NaCl, LiCl, Kr, and LiF. We find that the 
calculated bulk moduli by the hybrid functionals show better agreement with the experiments than 
the generalized gradient approximation (GGA) provides, whereas the calculated lattice constants 
by the hybrid functionals and GGA show comparable accuracy. The calculated energy band gaps 
and the valence-band widths for the ten prototype materials show substantial improvement using 
the hybrid functional compared with GGA. In particular, it is found that the band gaps of the ionic 
materials as well as the rare-gas solid are well reproduced by the LC-hybrid functional, whereas 
those of covalent materials are well described by the HSE functional. We also examine exchange 
effects due to short-range and long-range components of the Coulomb interaction and propose an 
optimum recipe to the short-range and long-range separation in treating the exchange energy. 



PACS numbers: 71.15.-m, 71.15.Mb, 71.20.Mq 
I. INTRODUCTION 

The local density approximation (LDA) (Ref. [l|) in 
density functional theory (DFT) (Ref. [3) has shown fan- 
tastic performance in understanding and even predict- 
ing material properties^ in spite of its relatively simple 
treatment of the exchange-correlation energy Exc[n] as 
a functional of the electron density n(r); e.g., for many 
materials, lattice and elastic constants are generally re- 
produced. The deviations from experimental values are 
within less than 1-2 % and several percent, respectively, 
in LDA. Yet LDA fails to describe some of properties 
including ground-state magnetic orderings even for bulk 
ironi and for some transition-metal oxides^ It also tends 
to overestimate the bonding strength, leading to an ab- 
solute error of molecular atomization energies^ 

Some of limitations of LDA are remedied by the gen- 
eralized gradient approximation (GGA) in which the 
exchange-correlation energy is expressed in terms of not 
only the electron density but also its gradient. The 
molecular atomization energies are calculated with the 
error of several tenths of eVj£ and the ground state of 
the bulk iron is correctly predicted to be a ferromag- 
netic body-center phased The prevailed functional form 
of GGA (PBE) (Ref. H) generally provides better accu- 
racy for structural properties of variety of solids and ac- 
tivation energies in chemical reactions than LDA does. 

The local (LDA) and semilocal (GGA) approximations 
are still insufficient to describe some of important prop- 
erties, however. The ground states of strongly correlated 



materials are incorrectly predicted and the energy band 
gaps of most semiconductors and insulators are substan- 
tially underestimated. The meta-GGA scheme^— ex- 
tending the exchange-correlation functionals with includ- 
ing the kinetic energy density further improves the LDA 
and GGA results for molecular systems^ but not succeed 
to remedy above failures in condensed matters. 

The failure of LDA and GGA is occasionally discussed 
in terms of the self- interaction error (SIE) i 13 i 14 An elec- 
tron is under the electrostatic potential due to other elec- 
trons. Yet the expression of the electrostatic potential 
in the (semi)local approximations includes the spurious 
interaction with the electron itself. When we consider 
the Hartree-Fock (HF) exchange potential with Kohn- 
Sham orbitals, this spurious self-interaction is cancelled 
by a term in the exchange potential. In the (semi)local 
expression of the exchange potential, however, this can- 
cellation is incomplete so that each electron is affected 
by the self-interaction. This SIE causes derealization of 
the electron, predicting the incorrect fractional-charged 
ground state of, e.g., with large nucleus separation^ 

The SIE affects the band gaps substantially. The 
band gap AE g is formally defined as the ionization en- 
ergy subtracted by the electron affinity so that Ai^g = 
E(N + 1) + E(N - 1) - 2E(N) where E(N) is the to- 
tal energy of the iV-electron system. In DFT with the 
exact exchange-correlation energy, the band gap is ex- 
pressed as the difference between the highest occupied 
Kohn-Sham level ejy+i(iV+l) of the (iV+l)-electron sys- 
tem and its counterpart of the iV-electron system sn(N): 
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i.e., AE g = e N+1 (N + 1) - e N (N)^tM When we intro- 
duce a fractional electron system with N + f electrons 
as a mixed state of real integer-electron systems, then 
the total energy E(N + /) becomes linear for < / < 1 
and shows discontinuity at the integer value N for finite- 
gap systems. Using Janak theorem^ which relates the 
Kohn-Sham level to the derivative of the total energy as 
e N+l (N + f) = dE(N + f)/df, the linearity of E(N + f) 
leads to the constant en + i(N + /) as a function of /. 
In the (semi)local approximations, however, the Kohn- 
Sham level £N+i{N+f) [£n(N— /)] increases (decreases) 
with increasing / due to the self-interaction, leading to 
the concave shape of E{N + /). This may cause an un- 
derestimate of the energy gap . 20 ' 21 

The HF approximation (HFA) is free from the self- 
interaction. Yet the calculated band gaps in HFA are 
substantially overestimated due to the lack of the corre- 
lation energy. An approach called the optimized effective 
potential 2 ^ which is incorporated in DFT (Refs. [23M25h 
intending to remedy the issue is still in an immature stage 
in a view of applications to polyatomic systems. 

Hence the hybrid functionals combining LDA or GGA 
with HFA may be effective to break the limitation of 
the semilocal approximations. The hybrid approach 
has begun in empirical ways: The HF-exchange energy 
was mixed with the LDA exchange-correlation energy 
in the half and half way 2 ^ and then three mixing pa- 
rameters were introduced^ 7 - to mix the LDA, GGA, and 
HFA energies; the latter scheme is called B3LYP and 
has widely been used to clarify thermochemical proper- 
ties of molecules i2£ A rationale for the hybrid functional 
is provided 2 ^ in the light of the adiabatic-connection 
theorem^ 

Exc N = f dX E xc .x , (1) 
Jo 

where 

£xc,a = (*A|K e |* A ) -~f*rj d^r'^^l (2) 

is the energy of the exchange and correlation in a 
system, where the electron-electron interaction V ee = 
(e 2 /2) — r j\ is reduced by the factor A but the 

external potential v\(r) is added to reproduce the elec- 
tron density n(r) of the real system (A = 1). Here is 
the ground-state many-body wave function. By assum- 
ing that -Exc. a is the 4th polynomial of A with particular 
asymptotic forms for A = and A = 1, Perdew, Ernz- 
erhof, and Becke have proposed a parameter-free hybrid 
functional called PBEO 2 ^ in which the HF and PBE ex- 
change energies are mixed with the ratio of 1:3. Its ap- 
plicability has been examined for molecular systems i21 

Screening of the Coulomb potential is effective in poly- 
atomic systems. Hence it may be effective to apply 
the non-local HF-cxchangc operator only to the short- 
range part of the Coulomb potential] 32 ' 33 This is con- 
veniently done by introducing the error function split- 



ting the Coulomb potential into short-range and long- 
range components. The hybrid functional which is con- 
structed in this way from the PBEO functional is pro- 
posed by Heyd, Sucseria, and Ernzerhof (HSE)^ 3 - This 
treatment reduces computational cost substantially and 
opens a possibility to apply the hybrid functionals to 
condensed matters. The structural properties as well 
as the band gaps of several solids have been calculated 
and significant improvements on semi-local functionals 
are achieved 

On the other hand, effects of the exchange interaction 
for the long-range component of the Coulomb potential 
are certainly importan t 39,40 in a view of reducing SIE. 
Hirao and his collaborators have proposed a long-range 
corrected (LC) functional in which the long-range com- 
ponent is treated by the HF exchange energy and the 
short-range component is by the LDA exchanged They 
have applied to various molecular systems and obtained 
relatively successful results^ 2 - - — Further application of 
the LC functional to molecular systems and its compari- 
son with other functionals have been done and the appli- 
cability of the LC functional have been recognized] 46 ' 47 
The LC functional has also been applied to structural 
properties and band gaps of several condensed matters 
and the results are compared with those obtained from 
other functionals4£ 

At the present stage, several hybrid functionals have 
been implemented in different packages and the assess- 
ment of the validity of each functional has been done 
mainly to molecular systems. Applicability of the PBEO 
and HSE functionals to condensed matters have been 
examined. 34 ' 35 ' 38 The structural and electronic proper- 
ties such as lattice constants, bulk moduli, and also the 
band gaps of condensed matters are obtained only after 
careful examinations of various calculation parameters. 
Obviously, numerical precision should not be neglected. 
In order to assess the validity of each hybrid functional, it 
is thus imperative to perform the computation in a single 
reliable calculation scheme. Furthermore, the split of the 
Coulomb potential into the short-range and long-range 
parts requires another parameter uj being the exponent 
of the error function. The w dependence of the results 
should certainly be examined for better understanding 
and further improvements on the hybrid functionals. 

The aim of the present paper is to implement sev- 
eral important hybrid functionals in the well-established 
plane- wave-basis total-energy band-structure calculation 
code and examine validity and limitation of each func- 
tional. A plane- wave code we adopt in this work is Tokyo 
Ab initio Program Package (TAPP)42r— We have cal- 
culated lattice constants, bulk moduli, band gaps, and 
band widths of various semiconductors and insulators 
with the PBEO, HSE, and LC hybrid functionals as well 
as (semi)local GGA functional. Comparison of the ob- 
tained results unequivocally elucidates the validity and 
the limitation of the hybrid functionals. 

In section [Til wc briefly describe each of the hy- 
brid functionals used in the present paper. Section IIIII 
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presents details of our computational scheme. The cal- 
culated results are shown in section ITVl and our finding 
is summarized in section Ivl 



II. HYBRID EXCHANGE-CORRELATION 
FUNCTION ALS 

In this section, we briefly describe the three hybrid 
functionals, PBEO, HSE, and LC, which we examine their 
applicability in this paper. 



A. PBEO functional 

Pcrdcw, Burke, and Ernzcrhof have proposed^ a poly- 
nomial form for -Exc, a as 
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where E% F and E^ FT are the exchange energies obtained 
by HFA and a certain (semi)local approximation in DFT, 
respectively. Here £ , xc'a=i * s the energy of the exchange 
and correlation defined as Eq. © and obtained by the 
(semi)local approximation in DFT. This formula infers 
that the (semi)local approximation in DFT is a good ap- 
proximation to Exc.x for A = 1. When A = 0, this 
formula equals to E^ F since i^xc^ = E^ FT for A = 0. 
From Eqs. ([TJ and we obtain 



Exc — E 



DFT 
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(4) 



Relying on the fourth-order Mollcr-Plesset perturbation 
theory applied for molecular systems, it is argued that 
n = 4 is the best choice.— Using the PBE functional^ as 
the approximation in DFT in Eq. (g]), the PBEO hybrid 
functional is given by 



Exc — E 



PBE 

xc 



\ { EW-E™), 



(5) 



leading to the mixing of 25 % HF exchange and 75 % 
PBE exchange. 



B. HSE functional 

Heyd, Sucscria, and Ernzcrhof have proposed 3 ^ a dif- 
ferent hybrid functional in which the long-range part of 
the HF-exchange energy is treated by the semilocal ap- 
proximation in DFT and the short-range part is calcu- 
lated exactly. The actual procedure is conveniently done 
by splitting the Coulomb potential as 



1 crfc(wr) erf(wr) 
r 



r 



r 



(6) 



and applying the first term only, i.e., the screened 
Coulomb potential, to the HF-cxchangc energy. The sec- 
ond term to the exchange energy is calculated with GGA. 



Adopting the mixing ratio in PBEO, then the HSE hybrid 
functional becomes 



pHSE 
Exc 



pPBE 
Exc 
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E x ' is the Fock-type double integral with the screened 
Coulomb potential. There is some complexity to divide 
the PBE exchange energy E^ BE into the short-range part 
£,pbe,sr an( j i on g_ ran g e p ar t £ ,FBE,LR . The dividing 
procedure will be shown in subsection MI CI 

Another ambiguous factor is the parameter uj in Eq. (|6]) 
which defines the short-range and the long-range parts of 
the Coulomb potential. The optimum value of uj = O.f 5 
a^ 1 (as'- Bohr radius) is proposed by examining the cal- 
culated results for molecular systems . 33 i 53 ~ — Examining 
the uj dependence of the calculated results for condensed 
matters is one of our aims in this paper. 



C. LC functional 

The HSE functional partly removes SIE by incorpo- 
rating the HF-exchange energy in the PBE functional. 
Yet the cancellation of the Hartrcc potential and the ex- 
change potential is absent in the long-range part. This 
may cause erroneous description of, e.g., the Rydbcrg 
states in isolated polyatomic systems or properties of 
charge-transfer systems. To remedy this point, appli- 
cation of the long-range part of the Coulomb potential 
to the HF-exchange energy is necessary. The long-range 
corrected (LC) functional has been proposed based on 
this viewpoint^ being expressed as 



pDFT 
Exc 



/ j-iHF.LR 
[Ex 



jtiDFT.LRn 
Ex )> 



(8) 



HF LR 

with E x ' being the Fock-type double integral with 
the long-range part of the Coulomb potential [the second 
term of Eq. JBJ]. For the DFT part, several approxima- 
tions including LDA^^, PBFJMSdL and other GGA 
(Ref. lill ) or meta-GGA (Ref. l4ot ) forms are adopted in 
Eq. and their validities are examined. It is argued 
that PBE combined in the LC hybrid scheme provides 
good accuracy for molecular properties^ As in the HSE 
scheme, there is some complexity to extract the long- 
range component of the exchange energy of the DFT 

, j-iDFT.LR 

part, E x 

The parameter uj in Eq. ^ affects the results sub- 
stantially. By examining the results for molecular sys- 
tems, the values ranging uj = 0.25-0.5 are argued to be 
optimum! 41 ' 46 By applying LDA in the LC-hybrid scheme 
to structural properties of solids, the value ui = 0.5 is 
found to produce reasonable results^ It is of interest to 
investigate the appropriate value of uj in the application 
of PBE in the LC-hybrid scheme to structural properties 
and band gaps of condensed matters. 
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III. COMPUTATIONAL DETAILS 

In this section, we describe our implementation of 
the hybrid functionals in the plane-wave-basis-set total- 
energy band-structure calculation code, TAPP42r— Nu- 
clei and core electrons are simulated by either norm- 
conserving^ or ultrasoftS pseudopotentials in the TAPP 
code. Non-locality of the HF-exchange potential gener- 
ally increases computational cost tremendously in the ap- 
plication to condensed matters. We have circumvented 
this problem using the Fast-Fourier transform (FFT), as 
explained below. The long-range nature of the Coulomb 
potential leads to singularity of its Fourier transform at 
the origin. This causes difficulty in numerical integration 
over the Brillouin zone (BZ) to obtain the HF-exchange 
energy in the PBEO and LC schemes. We here adopt a 
simple truncation scheme to overcome the problem. Fi- 
nally we explain how to divide the exchange energy in 
the PBE functional to the short-range and long-range 
components. 



A. Calculation of E% F by FFT 

The HF-exchange energy E^ F in condensed matters 
(crystal) are written as 



Ex_ F — n ^ ] / ^ /nk/n'k' Jnkn'k' ; 



(9) 



nk n'k' 



where f n k is the occupation number, and J n kn'k' is given 
by 



dirir ^" k(r) ^ k,(rV " k(rV "' k ' (r) . (10) 



Here </> n k( r ) is the Bloch-state orbital with the band 
index n and the wavevector k. The orbital <?W(r) is 
obtained by solving sclfconsistcntly the Euler equation 
(Kohn-Sham equation) in which the exchange-correlation 
potential is given by the functional derivative of the hy- 
brid exchange-correlation energy. The sum over k and n 
are taken for the occupied states. 

Our algorithm to compute the integral J„kn'k' in the 
plane-wave-basis code is as follows: Eq. (|10|) is written 
as 
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with u„k(r) = exp(— ik • r)0„k(r). In the plane-wave- 
basis-set scheme, the reciprocal lattice vectors G are 
used to represent the periodic wave function as u„k(r) = 
^ G e* G r M„k(G). We then obtain 
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When we compute Eq. (|T0|) directly, the calculation costs 
of J„ wk , and £HF are (jy3 ) and 0(iV b 2 and ) x 0(N£) x 
0(N G ), respectively, where Ng is the total number of 
the reciprocal vectors, Ak the total number of sampling 
k points in BZ, -/Vband the total number of the occupied 
bands. The number of Ng is much bigger than either 
^Vband or TVk- Hence the order N G is computationally 
demanding. 

We reduce this computational cost by using FFT. We 
first define the overlap density between states {nk) and 
(n'k') as 



"-nkn'k'(r) = w*k(r)M«'k'(r). 
Using this quantity, we obtain 

/^nkn'k' M^nkn'k' (r) 



(13) 
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We note that, since u„k(r) has unit-cell periodicity, the 
overlap density also has the same periodicity. Using the 
Fourier transformation of the overlap density n„k n 'k'(r), 
Eq. (fT4|) is rewritten as 



J, 



4- 



kn'k' 



^ Ik - k' 
G 1 



GP 



|n„ k „'k'(G)| 2 . (15) 



Since the n„kn'k'(G) can be calculated outside the sum- 
mation loop in Eq. (fT5|). the total calculation cost for 
J n kn'k' is 0(Ng + Ng^Ng)- The computational cost 
of E* F becomes O (A^ and ) x O (N* ) x O (N G + N G In N G ) . 
With considering the scaling as Ng oc N&tom, A^and oc 
Aatom, and N k oc l/N atom with iV a tom the number of 
atoms in the unitcell, the computational cost above is 
proportional to 7V a tom + A at om hi A a tom- 



B. Treatment of divergence in Coulomb interaction 

The Fourier transform v{q) of the Coulomb potential 
v{r) diverges at the long-wave length limit q — > 0. Calcu- 
lations of electrostatic energies thus require careful treat- 
ment and the well-known Ewald summation is the typical 
example. In calculating non-local exchange energies in 
condensed matters, more careful treatment is necessary 
We need to perform the BZ integration in evaluation of 
the exchange energy in Eq. ([9]). The integration is usu- 
ally performed by the summation with weighting factors 
of the integrand at finite discrete fe points, and thus en- 
counter a difficulty to evaluate the integral accurately 
by picking up the singular behavior of the Fourier trans- 
form of the Coulomb potential. There are several ways 
to overcome the difficulty. One is the auxiliary-function 
approach: An auxiliary function which has the same sin- 
gular behavior but is integrable is subtracted so that the 
summation can be done properly and the remaining term 
is obtained analytically— ~— An alternative way which 
we adopt in the present paper is simpler: We make the 
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Coulomb potential truncated at R c and then examine the 
convergence by numerically increasing R c £^ Namely, we 
replace Coulomb potential with a truncated potential, 



v(r) 



1F[ if \ r \ ^ R c 
otherwise 



(16) 



with a cutoff radius R c . What we need is of course con- 
verged quantities with R c — > oo. The truncated potential 
produces its non-divergent Fourier transform, 



v(q) 



4tt 



[1 - cos(|q|i? c )] 



(17) 



Convergence of required quantities with respect to R c 
is combined with the number of sampling k points in the 
BZ integration. The number of sampling fc points 
required should increase with increasing R c . We need to 
know the converged values with increasing both R c and 
Nk- The process to check this convergence in the two- 
parameter space can be done conveniently by introducing 
a certain relation between R c and Nk. We set up a rela- 
tion, Nk = 47ri?;?/(3f2 c ), where f2 c is the unit-cell volume, 
and examine the convergence of the exchange energy by 
increasing Nk and equivalently R c . We have found that 
the values of R c which are 4-6 times the dimension of the 
primitive unit cell are enough to assure the converged ex- 
change energies in ten materials calculated in the present 
paper. 



C. Calculations of i^ BE > SR and £ X BBE ' LR 

We next describe how to obtain the short-range and 
long-range parts in the PBE-exchange energy, following 
Ref. |62|. The original expression for the PBE-exchange 
energy is given by, 



E 



PBE = / drer if [n(r)]n(r)F x PBE (s) , (18) 



where e x [n] is the exchange energy of the homogeneous 
electron gas with the electron density n and F PBE (s) 
is an enhanced factor due to a density gradient s = 
V?i|/(2fcpn) with the Fermi wavenumber fcp = (i^n) 1 / 3 . 
The enhanced factor is written in an integral form 



PBE/ 
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(s) = --l dyyJ™»(8,y), 



(19) 



where y = kpr is a dimensionless quantity and J PBE (s, y) 
describes an exchange-hole density at the distance r. 
Heyde, Scuseria, and Ernzerhof^ proposed an expres- 
sion for the short-range PBE-exchange functional, where 
the original Coulomb interaction 1/r is modified to a 
screened form erfc(ajr)/r. This modification leads to an 



enhanced factor somewhat different from the original ex- 
pression of Eq. (fTi?]) as 



F PBE ' SR (s, W ) = -^ j( ^J PBE ( S , y )erfc^J.(20) 

The short-range exchange energy is given with this en- 
hanced factor as 

£PBE,SR (a;) = J dre unif [n(r)]n(r)F PBE,SR (S)a;) , (21) 

The long-range PBE-exchange term is defined by sub- 
tracting the short-range part in Eq. (|2~T|) from the original 
one in Eq. (fT5|) as 



E^ lh (uj) = - E^-^iuj). (22) 

Implementation details for these calculations can also be 
found in Ref. I&a. 



D. Calculation conditions 

We generate norm-conserving pseudopotential to sim- 
ulate nuclei and core electrons, following a recipe by 
Troullier and Matins^ The core radius r c is an essential 
parameter to determine transferability of the generated 
pseudopotential. We have examined r c dependence of the 
calculated structural properties of benchmark materials 
and adopted the pscudopotentials generated with the fol- 
lowing core radii in this paper: 0.85 A for Si 3s, and 1.16 
A for Si 3p, 1.06 A for Ge 4s and 4p, 1.06 A for Ga 4s 
and 4p, and 1.48 A for Ga Ad, 0.64 A for N 2s and 2p, 
0.85 A for C 2s and 2p, 0.79 A for O 2s and 2p, 1.38 A 
for Na 2s and 2p, 1.16 A for CI 3s and 3p, 0.95 A for Li 
2s, 0.64 A for F 2s and 2p, 1.59 A for Ba 5s, 5p, and 5d, 
1.32 A for Mg 2s and 2p, and 1.38 A for Ti 3d and 4s, 
and 1.43A for Ti 4p, 1.48 A for Kr 4s and 1.37 A for Kr 
4p. 

The pscudopotentials are generated by the (semi) local 
approximations in DFT. This means that the HF- 
exchange energy between core and valence states are ne- 
glected. Yet we have found that the calculated energy 
bands obtained with the pscudopotentials generated in 
LDA and GGA, are essentially identical to each other, 
implying that the treatment of the exchange-correlation 
energy in generating pscudopotentials has minor effects. 

The partial core correction^ is not included in our cal- 
culations. This is partly because magnetic properties are 
not considered in the present work. However, to assure 
the accuracy of structural properties, we regard some of 
core orbitals as valence orbitals and include them ex- 
plicitly in the pseudopotential generation: Such orbitals 
included as valence states are 2s and 2p orbitals of Na 
and 2s and 2p orbitals of Mg. 

Appropriate choice of cutoff energies E cut in the plane- 
wave-basis set, which is related to hardness of the 
adopted norm-conserving pscudopotentials, is a princi- 
pal ingredient to assure the accuracy of the results. We 
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(a) PBE 



(b)HF 



(c) PBEO 




XK r L KWX 
(d) HSE (co = 0.1) 



XK r L KWX 
(e) HSE (co = 0.2) 



XK r L KWX 
(f) HSE (co = 0.3) 




XK T L KWX 
(g)LC (co = 0.1) 



XK r L KWX 
(h) LC (co = 0.2) 



XK r L KWX 
(i) LC (co = 0.3) 




KWX 



FIG. 1: (Color online) Calculated energy bands of Si with PBE (a), HF (b), PBEO (c), HSE (d)-(f), and LC (g)-(i) exchange- 
correlation functionals. The origin of the energy is set at the valence band top. 



have examined convergence of structural properties and 
band gaps with respect to E cut and reached the following 
well converged values with E cut for each material: 25 Ryd 
for Si and Ge, 36 Ryd for Kr, 64 Ryd for BaTi0 3 , and 
100 Ryd for diamond, GaN, MgO, NaCl, LiCl, and LiF. 
The remaining important ingredient to assure our assess- 
ment of each hybrid functional is the sampling k points 
for the BZ integration. We have adopted the scheme by 
Monkhorst and Pack in which BZ is divided by equally 
spaced mesh. After careful examination, we have found 
that 4x4x4 sampling k points are enough to assure the 
accuracy of the total energies and energy bands in the 
ten materials. The results are confirmed by repeating 
the calculations with 6x6x6 sampling k points. 



IV. RESULTS AND DISCUSSION 



We have performed total-energy electronic-structure 
calculations using PBE, HF, PBEO, HSE, LC exchange- 
correlation functionals for ten materials, including co- 
valent semiconductors, ionic insulators, dielectric com- 
pounds, and a rare-gas solid: i.e., Si, Gc, GaN, diamond, 
MgO, NaCl, LiCl, LiF, BaTi0 3 , and Kr. The calculated 
results elucidate capability and limitation of each func- 
tional in discussing electronic and structural properties 
of these prototype materials. We first present calculated 
electron states of these materials and compare them with 
experimental results in the subsection IIV Al Then we 
present the results of structural optimization in the sub- 
section IIVB1 
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(a) PBE (b) HF ( C ) PBEO 




r m x r r tmxt r rMxr r 



(d)HSE (co = 0.1) (e) HSE (to = 0.2) (f) HSE (co = 0.3) 




r m x r r r m x r r rivixr r 

(g) LC (tu = 0.1) (h) LC (co = 0.2 ) (i) LC (co = 0.3 ) 




rivixr r rMxr r rMxr r 



FIG. 2: (Color online) Calculated band energy bands of BaTi0 3 with PBE (a), HF (b), PBEO (c), HSE (d)-(f), and LC (g)-(i) 
exchange-correlation functionals. The origin of the energy is set at the valence band top. Note that the vertical scale in (b) is 
different from those in other panels. 



A. Energy bands and gaps 

Figure [T] shows calculated band structures of Si with 
five exchange-correlation functionals. For the HSE and 
LC functionals, the band structures with different choices 
of w [0.1, 0.2, and 0.3 (a^ 1 ) ] in Eq. © are shown. The 
overall features of the band structures obtained by the 
five functionals are similar to each other. Yet the band- 
widths and the fundamental gaps arc different quantita- 
tively. When the ratio of the HF exchange to the total- 
exchange functional is large, the resulting bandwidth and 
gap become large; these quantities become larger in the 
order of PBE, PBEO, HSE, LC, and HF. Notice that 
the band gaps obtained by the LC functional are always 
larger than those by the HSE functional, indicating that 
the correction to the long-range part of the exchange po- 
tential tends to make the band gap large. 

We also show the calculated energy bands of the dielec- 
tric compound BaTiOa and the ionic insulator NaCl in 
Figs. f2] and [31 respectively. We find the general tendency 



similar to that in Si: i.e., the overall features of the energy 
bands arc insensitive to the difference in the functionals; 
the LC functional provides larger energy gaps compared 
with the HSE functional. 

TABLE HI and HT1 summarize the calculated bandwidths 
and band gaps for the ten materials. For assessment of 
validity of each functional, it is convenient to introduce 
two quantities, the mean relative error (MRE) which is 
the mean of the calculated value minus the experimen- 
tal value over the ten materials, and the mean absolute 
relative error (MARE) which is the mean of the absolute 
value of the difference between the calculated and the 
experimental values over the ten materials. The MRE is 
a measure of under- or over-estimates of the experimen- 
tal values since each functional predicts either smaller 
or larger values than the experimental values for most 
of the ten materials. On the other hand, MARE is a 
measure of closeness between the calculated and experi- 
mental values. In the discussion below, we categorize the 
ten materials into Group I and Group II: Group I, i.e., 
diamond, GaN, BaTiC^, Si and Ge has covalent charac- 
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(a) PBE 



(b)HF 



(c) PBEO 




r k 

(d)HSE (o)=0.1 ) 



r k 

(e) HSE (00=0.2 ) 



K 



(f) HSE (co=0.3 ) 




FIG. 3: (Color online) Calculated energy bands of NaCl with PBE (a), HF (b), PBEO (c), HSE (d)-(f), and LC (g)-(i) 
exchange-correlation functionals. The origin of the energy is set at the valence band top. 



ters in which the band gaps are less than 7 eV; Group II, 
i.e., LiF, Kr, NaCl, LiCl and MgO, consists of ionic solids 
and a rare-gas solid in which the experimental band gaps 
are larger than 7 eV. 

The calculated band gap by the PBE functional for 
each material is substantially smaller than the corre- 
sponding experimental value, as is reported in literature. 
The calculated MRE for the ten materials is —42.5 %. 
On the other hand, HFA largely overestimates the band 
gaps for all the ten materials: MRE is 179.7 %. The 
PBEO functional provides better values. The calculated 
MRE by the PBEO for the materials in Group I is 25.4 
%, whereas it gives —14.0 % for the materials in Group 
II. 

The HSE and LC functional also provide better agree- 
ment of the band gap with the experiments than the PBE 
does. Further when we choose optimum values of w, the 
HSE and LC results show nice agreement with the exper- 
imental values. In the HSE functional, a general trend is 
the decrease in the band gap with increasing u>. We have 
found that the value us = 0.1 a^ 1 produces the band gaps 



close to the experimental values within MRE = 0.9 % for 
the materials in Group I. It produces worse values for 
the materials in Group II. The calculated MRE is —19.0 
% for the materials in Group II, but still HSE being the 
better approximation than PBE. 

The LC functional provides good agreement with the 
experimental values for the materials in Group II: When 
we choose the optimum value of w = Q.2a^ B 1 1 the cal- 
culated MRE is nicely small, 4.4 % for the materials in 
Group II. Yet the LC provides worse values for the mate- 
rials in Group I, showing its limitation as an universally 
valid approximation. 

Figure U is a summary of our calculated band gaps by 
the PBE, HSE and LC functionals. For HSE and LC, we 
show the calculated results with the optimum lu values: 
lo = 0.1 cig 1 for HSE and w = 0.2 o^ 1 for LC. It is clearly 
shown that the calculated band gaps by hybrid function- 
als, HSE and LC, are in better agreement with the exper- 
imental values than the PBE (GGA) approximation, in- 
dicating the promising possibility of the hybrid function- 
als. The degree of the agreement is close to the Green's- 
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TABLE I: Band gaps e gap obtained from PBE, HF, PBEO, HSE, and LC calculations. The uj is a parameter which separates 
the long-range and short-range parts of the Coulomb interaction (see text) . Experimental values are taken from Ref. HH for Ge, 
Si, and C, Ref. H for /3-GaN, Ref. [7(3 for BaTiOs , Ref. [7l| for MgO, Ref. El for NaCl, Ref. [H for Kr, and Ref. [fl for LiCl, and 
LiF. The calculated mean relative error (MRE) and the mean absolute relative error (MARE) with respect to experimental 
values are also shown in percent. Group I consists of materials having the experimental gap less than 7 eV, while Group II of 
materials with the gap more than 7 eV (see text). 
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56.3 


14.0 


19.0 


23.9 


27.3 


29.8 


15.6 


5.6 


18.9 
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56.3 
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TABLE II: Valence bandwidth W obtained from PBE, HF, PBEO, HSE, and LC calculations. The uj is a parameter which 
separates the long-range and short-range parts of the Coulomb interaction (see text). Experimental values are taken from 
Refs. EKzi for Ge, Refs. \EMM for Si, Refs. EMM!! for C. 
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7.44 
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1.67 
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2.38 
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2.82 
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3.39 


3.16 


3.16 
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3.59 


3.91 


4.17 




Kr 


1.50 


1.93 


1.61 


1.59 


1.55 


1.52 


1.57 


1.74 


1.85 




LiF 


2.83 


3.66 


3.09 


2.96 


2.92 


2.87 


2.94 


3.09 


3.27 





function-based GW approximation in which there are 
several ambiguities in theoretical treatments £5r— How- 
ever, this figure also shows the limitation; HSE is reason- 
ably good for only the Group I materials, whereas LC is 
good for only the Group II materials. 

Our finding here is that HSE is a good approxima- 
tion for relatively small-gap materials and that LC is a 
good approximation for the relatively large-gap materi- 
als. Screening of the Coulomb interaction depends on 
the materials. Hence the best treatment of the short- 



range and long-range parts of the Coulomb interactions 
should change material by material. What we have found 
here is natural in the sense. However, another point we 
have shown here is that the appropriate choice of uj with 
each exchange-correlation functional provides reasonable 
agreement with the calculated band gaps for rather wide 
range of materials; HSE with uj ~ 0.1 cig 1 for the ma- 
terials with their band gaps smaller than 7 eV and LC 
with uj ~ 0.2 a^ 1 for the materials with their band gaps 
larger than 7 eV. This information certainly provides a 
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practical recipe to obtain reliable band gaps for various 
materials. 

The calculated bandwidths by the PBE functional are 
in good agreements with the experimental values avail- 
able (TABLE [II]). Other hybrid functional also pro- 
vide reasonable agreement, in particular, with appro- 
priate choices of the 10 parameter for the HSE and LC 
functionals. On the other hand, HFA shows substantial 
overestimate for each material. 



B. Structural properties 

We next examine performance of the hybrid function- 
als in describing structural properties. TABLE IIIII and 
IIVI show calculated lattice constants ao and bulk moduli 
-Bo, respectively, of the ten materials. The ao and Bq 
values are determined by fitting parameters in the Mur- 
naghan equation of state to the calculated total energy 
as a function of the volume. In MARE and MRE pre- 
sented in TABLE IIIII the calculated value for solid Kr- 
is not included. As shown in this TABLE, the calcu- 
lated lattice constant of Kr is substantially larger than 
the experimental value: It is overestimated by 1.1-9.2 % 
in GGA and in the hybrid approximations, and by 1.2 % 
in HFA. This is because van der Waals interaction which 
is unable to be treated in the approximations examined 
in this paper plays an essential role in solid Kr. We thus 
exclude the value for Kr from the statistical assessment. 
The results obtained by PBE reasonably agree with the 
experimental values (MRE = 0.68 % and MARE = 1.20 
% for a and MRE = -7.1 % and MARE = 8.3 % for 
Bq). The accuracy of HFA is slightly inferior to that of 
PBE (MRE = -1.10 % and MARE = 1.37 % for a and 
MRE = 13.7 % and MARE = 13.7 % for B ). The PBEO 
functional shows better accuracy with the MRE = —0.49 
% and MARE = 1.10 % for o , and MRE = -1.7 % and 
MARE = 5.1 % for B . 

The calculated lattice constants using the HSE and LC 
functionals are insensitive to the choice of uj: The differ- 
ence obtained from different uj values is within 1 %. The 
difference in the calculated bulk modulus from different 
uj values are not so small as in the lattice constant partly 
because the fitting by the Murnaghan equation is incom- 
plete. For the HSE functional, the value u = 0.1 Og 1 
produces best agreement with the experiments: MRE = 
-0.47 % and MARE = 1.40 % for a , and MRE = 1.7 
% and MARE = 3.6 % for B . For the LC functional, 
we have found that uj = 0.2 o^ 1 produces best results: 
MRE = 0.10 % and MARE = 1.21 % for a , and MRE 
= 2.0 % and MARE = 5.5 % for B . These optimum 
values for uj for the HSE and LC functionals are identical 
to the optimum values determined from the calculated 
MRE and MARE for band gaps, corroborating the ap- 
propriate choice of uj. 



V. CONCLUSION 

We have studied validity of hybrid exchange- 
correlation functionals in density functional theory 
by implementing three hybrid functionals in a well- 
established plane-wave-basis-set code named TAPP and 
by calculating structural properties and electron states 
of representative ten materials where the experimental 
energy gaps range from 0.67 cV to 14.20 cV. The three 
hybrid exchange-correlation functionals examined in this 
paper are PBEO proposed by Perdew, Burke, and Ernz- 
crhoff, HSE proposed by Heyd, Scuseria, and Ernzerhoff, 
and LC originally proposed by Savin and by Hirao and his 
collaborators. For comparison, Results from the gener- 
alized gradient approximation, i.e., PBE, and HFA have 
been presented. The ten materials we have examined arc 
Ge, Si, GaN, BaTi0 3 , diamond, MgO, LiCl, NaCl, Kr, 
and LiF which are representatives of covalcnt, ionic, and 
rare-gas solids. 

We have found that the structural properties such as 
lattice constants are already well reproduced by the PBE 
functional and also by HFA and that the hybrid function- 
als show better agreement with the experimental values. 
We have determined appropriate values of uj in the sep- 
aration of the short-range and long-range parts in the 
Coulomb interaction: The optimum value is uj = 0.1 a^ 1 
for the HSE functional and u> = 0.2 a^ 1 for the LC func- 
tional. By choosing the appropriate value of u> in the 
HSE and LC functionals, we have achieved better agree- 
ment in the lattice constants and further substantial im- 
provement in description of elastic constants such as bulk 
moduli for the ten materials. 

Dramatic success of the hybrid functionals are ob- 
served in the calculated band gaps. We have found that 
the calculated band gaps by the LC functional for the 
wide band-gap materials satisfactorily agree with the ex- 
perimental values with mean relative error (MRE) of 3.0 
%, whereas the band gaps by the HSE functional for the 
small band-gap materials agree well with the experimen- 
tal values with MRE = 0.7 %. This good description of 
the band gaps is unprecedented in density functional the- 
ory where LDA and GGA produce the value of MRE ~ 
40-50 %, and is comparable with or better than what the 
GW approximation produces. The u value leading to the 
best agreement with the experiments is 0.1 a^ 1 for the 
HSE functional and 0.2 a^ 1 for the LC functional. These 
uj values are identical to the optimum values determined 
from the examination of the structural properties. The 
calculated valence-band widths by the hybrid functionals 
also agree satisfactorily with the experimental values. 

It is now established that the HSE and LC functionals 
with appropriate choice of the uj parameter are useful to 
describe structural and electronic properties of various 
materials. Rigorous justification of the choice of the form 
of the hybrid functionals along with a guiding principle 
of choice of uj would offer further developments in the 
first-principles calculations. 
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TABLE III: Lattice constants ao (A) obtained from PBE, HF, PBEO, HSE, and LC calculations. The uj parameter separates 
the short-range and long-range parts of the Coulomb interaction (see text). Experimental values are taken from Ref. [6^ for 
Ge, Si, /3-GaN, C, and MgO, Ref . [zlzl for BaTiOs, Ref. for NaCl, LiCl, and LiF, and Ref. HD for Kr. The calculated mean 
relative error (MRE) and the mean absolute relative error (MARE) with respect to experimental values are also shown in per 
cent. 
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MARE (%) 
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TABLE IV: Bulk moduli Bo (GPa) obtained from PBE, HF,PBE0, HSE, and LC calculations. Experimental values are taken 
from Ref. [H for Ge, Ref. for Si, /3-GaN, C, MgO, NaCl, LiCl, and LiF, Ref. for BaTiOs, and Ref. HD for Kr. The 

calculated mean relative error (MRE) and the mean absolute relative error (MARE) with respect to experimental values are 
also shown in per cent. 
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FIG. 4: (Color online) Calculated band gaps obtained from different exchange-correlation functionals: PBE (blank squares), 
HSE with w=0.1 ag 1 ( green dots), and LC with uj—Q.2 a B (purple squares), plotted against experimental band gaps. Group 
I consists of materials having the experimental gap less than 7 eV, while Group II is composed of materials with the gap more 
than 7 eV. 
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